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We propose a simple feedback-control scheme for adiabatic quantum computation with supercon- 
ducting flux qubits. The proposed method makes use of existing on-chip hardware to monitor the 
ground state curvature, which is then used to control the computation speed in order to maximise 
the success probability. We show that this scheme can provide a polynomial speed-up in perfor- 
mance and that it is possible to choose a suitable set of feedback-control parameters for an arbitrary 
problem Hamiltonian. 



I. INTRODUCTION 

Adiabatic quantum computing (AQC) is a promising 
paradigm of quantum computation. It was first proposed 
as a possible method of solving NP-complete problems in 
a more efficient manner than is possible classically [l| . It 
has since been shown to be equivalent to other univer- 
sal models of quantum computation 0, Q . There have 
also been some significant advances in the experimental 
implementation of AQC using superconducting flux 
qubits 0,[||. 

In AQC, the ground state of a Hamiltonian encoding 
the problem is reached by adiabatic evolution from a con- 
figuration with an easily-reachable ground state. The fi- 
nal ground state will encode the solution to the problem. 
This process can be described by a Hamiltonian of the 
form 



H [A(t)] =H P + X(t)H b 



(1) 



where H p describes the problem, iff, describes a strong 
applied bias (such that Hb 3> H p and [Ht,,H p ] ^ 0) 
and the bias strength, X(t), is slowly varied from 1 to 0. 
The instantaneous energy eigenvalues, Et (A), and eigen- 
states, |^(A)}, of ^ for an n-qubit system will be given 
by 



occur towards the end of the computation, when the lev- 
els become more densely packed as the relative strength 
of the perturbing Hf, is much less than H p , as seen in 
[g.llCj|. In the vicinity of this avoided crossing, the chance 
of excitation through Landau-Zener-Stueckelberg (LZS) 
[ill ] tunnelling is at its maximum. The chance of exci- 
tation via the LZS mechanism increases monotonically 
with dX/dt as well as A m j n . As A m j n is fixed by the 
choice of problem and bias Hamiltonians, the rate dX/dt 
must be chosen to be relatively slow in order to minimise 
the excitation probability. The optimum speed to cross 
the critical minimum gap in the presence of decoherence 
was derived in [T2]. This is the optimal speed for the 
entire computation, if it is performed at a constant rate. 
However, far away from the minimum gap, where the 
ground and first excited states are well separated, the 
computation could be run much faster than this optimal 
speed without significantly increasing the chance of exci- 
tation. This leads to the idea of designing time-optimal, 
non-linear interpolation paths, X(t), for AQC that take 
into consideration the critical points in the evolution. A 
number of approaches to realize this have been explored; 
13 propose a variational ansatz that allows an op- 
timal path to be determined, and [l5| suggests that de- 
phasing may pick out the optimal path under certain 
conditions. 



H(X) \l (A)) = E t (A) \i (A)) 



(2) 



where < I < 2™ - 1 and E < E x < . . . < £ 2 »-i- 
In most AQC schemes, it is usually assumed that the 
interpolating parameter X(t) varies linearly with time and 
therefore the rate dX/dt = constant. 

The critical point in the evolution of an adia- 
batic quantum computation is the minimum gap be- 
tween the ground and first excited state, A m j n = 
min 0<A ( t )< 1 [Ei(X) — E (X)]. For an AQC process de- 
scribed by a total Hamiltonian of the form ([T]) , A m i n will 
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In this paper we propose a feedback-control scheme for 
the interpolation rate based on the result of weak con- 
tinuous measurements of the ground state curvature of a 
flux qubit AQC system. The aim is to allow the compu- 
tation process to move quickly when the ground state is 
well separated from the excited states, but then gradu- 
ally slow down as it approaches and traverses the critical 
minimum gap. After it has moved through the minimum 
gap, the rate of the computation can then be safely in- 
creased again. We show that this type of control scheme 
can provide a significant improvement in the performance 
of existing solid state implementations of AQC without 
the need for additional on-chip hardware. 
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II. SPECIFICATIONS FOR 
FEEDBACK-CONTROL OF AQC 

Feedback-control [HI ] is a very common control scheme 
that has been applied to a wide range of applications. Es- 
sentially, this approach makes use of the results of contin- 
uous measurement of a dynamical system to continuously 
adjust its control parameters to regulate its output. 

In terms of AQC, the main process variable that we 
would like to control is the probability of successfully 
finding the system in the ground state of the problem 
Hamiltonian at the end of the evolution, when TL(X = 



0) = Hp. This success probability is given by 
P=|(^(A = 0)|0(A = 0))| 2 



(3) 



The aim of a feedback-control scheme for AQC would be 
to increase P without the need to also increase the to- 
tal computation time T. The most easily-controlled pa- 
rameter in (TTJ) is the interpolation function X(t) and its 
corresponding rate of change dX/dt. With relatively pre- 
cise control of dX/dt throughout the evolution, it would 
be possible to move slowly through the critical minimum 
gap and then evolve quickly when the ground state is 
well separated from the rest of the spectrum. This type 
of scheme should in theory allow a given probability to be 
achieved with a shorter T when compared to the linear 
interpolation of A(t), and therefore produce an improve- 
ment in the scaling of P with T. 

The key requirement needed to realise this scheme is a 
means of detecting when the system is moving through 
the critical minimum gap. Crucially, the measurement 
must have minimal back-action on the evolution of the 
system. It also needs to be continuous and allow read- 
out on a timescale that is much shorter than the total 
computation of the AQC operation. 

A method of weak continuous measurement, that has 
been succesfully applied to systems of superconducting 
flux qubits, is the impedance measurement technique 
(IMJ), see @ for a detailed review of IMT or e.g. 
[1, [13, El] for applications of it. In the IMT, a driven 
high-quality LC tank circuit is inductively coupled to 
the flux qubit system. The coupling to the qubits will 
induce a frequency shift in the resonator which is pro- 
portional to the curvature of the qubit system's state. 
This frequency shift can be monitored by measuring the 
impedance of the tank through the phase angle between 
the driving current and the tank voltage. In AQC, the 
system should remain in the instantaneous ground state 
throughout the evolution and therefore the results of the 
IMT will be proportional to the ground state curvature, 
d 2 Eo(X) /dX 2 . The ground state curvature is a particu- 
larly appropriate property of the system to monitor in a 
feedback-control loop, as we know that it will increase to 
its global maximum at the critical point A m j n , as shown 
in the example of Fig. [TJ 



_= 

> 

a 

fct) 
W 




100 



50 



> 

u 



0.5 

Bias strength A.(f) 



FIG. 1. (Color online) Plot of the lowest two energy eigenval- 
ues, Eo(X) and -Ei(A) (black solid lines), and the curvature of 
the ground state, d 2 Eo(\) /d\ 2 (red dashed line), for a typical 
example of a 2-qubit adiabatic quantum computation. 



The back-action of the IMT measurement on the qubit 
system was studied in detail by Smirnov in [l9j . The 
conditions required to consider IMT as a weak continuous 
measurement of Rabi oscillations in a flux qubit were 
determined and it was found that the contribution to the 
decay rate of the qubit due to the back-action of the tank 
circuit is peaked when = lolc, where is the qubit 
Rabi frequency and u>lc is the resonant frequency of the 
tank circuit. In this case the relevant frequency scale of 
the qubit system is A m - ln /h and we therefore require 



> CJ LC + 1LC 



< U LC - JLC, 



(4) 



where jlc is the tank lincwidth, in order to minimise 
the decoherence caused by the IMT measurement and 
allow it to be classified as weak continuous measurement. 
Requiring the two frequency scales to be out of resonance 
is a relatively easy requirement to meet practically in 
comparison to requiring them to be in resonance. This is 
because it still does not require exact knowledge of A m i n , 
only an order of magnitude estimate. 

Therefore, IMT is a promising potential candidate 
for the type of weak continuous measurement required 
to practically implement a feedback-control scheme for 
AQC, where we take 



dX 
dt 



d 2 E (X) 
dX 2 



(5) 



The IMT also has the advantage that the same LC tank 
circuit can also be used for dispersive readout of the so- 
lution state, as in Q, which means that the proposed 
feedback-control scheme does not require any additional 
on-chip hardware. 

We envision two possible modes of operation for this 
scheme in experimental implementations: (i) Offline con- 
trol. Where one initial evolution of the system is per- 
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formed with the weak continuous measurement to pre- 
computc d 2 Eo(\)/d\ 2 . Then the actual AQC evolution 
is carried out using the resulting interpolation function as 
many times as necessary to collect the required statistics 
for the measurement of the solution state, (ii) Online 
control. In this case the weak continuous measurement 
is performed throughout all of the AQC evolutions and 
the rate is adjusted actively through the feedback-control 
loop. In this case, the effects of the back-action of the 
measurement on the system will be a more important 
consideration. 



by 



T = 



dX. 



(10) 



IV. A PROTOTYPICAL EXAMPLE OF 
FEEDBACK-CONTROLLED AQC 

We take the bias Hamiltonian Hb to be a strong trans- 
verse field acting equally on all n qubits with an easily 
constructed Hamiltonian of the form 



III. THEORETICAL MODEL FOR 
FEEDBACK-CONTROLLED AQC 

We assume that the inequality Q holds true and anal- 
yse the ideal case of a negligibly small back-action from 
the continuous IMT measurement on the qubit system. 
In order to simulate this system, we first need to ac- 
curately compute the ground state curvature as a func- 
tion of A(t). The Pechukas- Yukawa equations of motion, 
BUS US can be used to describe the dynamics of the 
energy eigenvalues of a quantum system with a Hamilto- 
nian of the form ([!]); 



dE t 
~dX 

dve 
dX 

dX 



E 



21/ 



Ck-\ 



(6) 



^2 liklk 3 



{E t - E k y (Ej - E k y 







(11) 



where the constant Z»l and a a (with a = x,y,z) are 
the usual Pauli matrices. The ground state |0&) of Hb 
is simply an equal superposition of all 2™ computational 
basis states, 



2™-l 



|o h ) = 2-/ 2 J2 li) 



(12) 



In the following calculations we take the Z = 1(W 2 . 

As a prototypical example of an AQC problem Hamil- 
tonian we use a H p with randomly weighted coupling in 
the z-dircction between all possible combinations of the 
n qubits; 



H, 



= E 



(13) 



where jt is the i bit in the binary representation of j 
and tj are the coupling strengths, which are selected from 
where v e = (£\ H b \£) and l tj = (E e - E 3 ) (£\ H b \ Then a Gaussian distribution with mean /i = and standard 
from ©, the ground state curvature will be given by 

2 E -21/qiI 2 



deviation a 



This form of randomly selected prob- 



dX 2 



(Ei - EoY 



(7) 



The evolution of the state of the system, \ip (A)), will 
be described by the Schrodinger equation 



H(X)\^(X)) 



(8) 



The (dX/dt) ter m in ([5J describes the feedback con- 
troller rule and will, in general, be of the form ([5]). The 
simplest type of feedback controller is a constant propor- 
tional gain controller and, as we wish the interpolation 
rate to decrease as the curvature increases near the crit- 
ical minimum gap, we therefore take 



dX 

~dt 



=^=ik 



d 2 E 
dX 2 



(9) 



where k is the gain constant. The total computation time 
T of a feedback-controlled AQC evolution will be given 



lem Hamiltonian represents a reasonable approximation 
as it can be used to encode any finite computational op- 
timisation problem with a proper choice of 6j. 

To investigate the efficacy of the proposed feedback- 
control scheme for AQC, we have solved equations (0 
and ([S]) numerically for a range of different gain constants 
k and system sizes n for an ensemble of random i/ p 's. We 
also solve © for a range of linear interpolations, where 
dX/dt = constant, so that we have a benchmark for com- 
parison to the results for feedback-controlled AQC. In 
both cases the system is initially prepared in the ground 
state of the total Hamiltonian ([1]), |0 (A = 1)). 

A comparison of the scaling of the success probability 
P with the total computation time T for a typical 2- 
qubit AQC operation with and without feedback-control 
is shown in Fig. In this case, the feedback-controlled 
interpolation causes P to approach unity more rapidly 
than the standard linear interpolation scheme allows. 
This improvement can be interpreted in one of two ways; 
as an increase in P at a given T, or conversely as a de- 
crease in the T required to achieve a given P. 
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FIG. 2. (Color online) Comparison of the scaling of the suc- 
cess probability P versus the total computation time T for lin- 
ear and feedback-controlled interpolation for a typical exam- 
ple of a 2-qubit adiabatic quantum computation. The dashed 
black line indicates the adiabatic time, T a( j = 2.775, for the 
computation. The inset shows the non-adiabatic oscillations 
that arise for the feedback-controlled case around the adia- 
batic limit. 



Figure [5] also shows that the success probability, P, for 
the feedback-controlled interpolation oscillates slightly in 
the vicinity of the adiabatic time p2j . 



'ad 



max 

1<j<2"-1 



max 

0<A<1 



|Q(A)|ffi|0(A))| 



(14) 



These non-adiabatic oscillations most likely arise from 
the system's non-linear trajectory through phase space, 
and despite them the feedback-control scheme still offers 
a significant improvement over linear interpolation. In 
the adiabatic limit, when T 3> Tad, feedback-controlled 
interpolation will still offer some improvement, but as P 
will already be arbitrarily close to 1, this improvement 
will be negligible. The regime where the feedback-control 
method excels is the approximately adiabatic case, where 
T ~ T a d, as can be seen in Fig. [2J In practical AQC 
devices we would like to achieve a relatively high suc- 
cess probability P in the shortest time T possible and 
it is the approximately adiabatic regime that offers the 
largest potential gains in this trade-off. This is because 
a small increase in T can significantly increase P and the 
proposed feedback-control scheme improves this situation 
even further. 

A key issue for our proposed feedback-control approach 
is how well its performance scales with the system size. 
The dependence of the average computation time on the 
system size for linear and feedback-controlled interpola- 
tion is shown in Fig. [3] For linear interpolation we have 
T = 0(n 3 3 ), while for the feedback-controlled interpola- 
tion we find that T = 0[n 1A ) . Therefore, in this case the 
proposed feedback-control method provides a polynomial 



FIG. 3. (Color online) Comparison of the average compu- 
tation time T required to achieve a success probability of 
P = 90% for linear and feedback-controlled interpolation 
as a function of the number n of qubits. The polynomial 
lines of best fit give T = 0(n 3 ' 3 ) for linear interpolation and 
T — 0(n 1 ' 4 ) for the feedback-controlled case; an approxi- 
mately quadratic speed-up. The results are averaged over 
10,000 randomly selected problem Hamiltonians of the form 
(|13|l and the error bars denote the standard deviation of the 
distributions of T. 



because as n increases the system's spectra become more 
densely packed and A min decreases (as shown in [23j|). 
This means that the speed at which A m j n is traversed 
becomes an even more important factor in determining 
whether the system is excited out of the ground state. 
Therefore, controlling the speed in proportion to the gap 
becomes an even more effective strategy. We note that 
it is impossible to infer exactly how the maximum T will 
scale in the limit of large n from these results though. 



Another consideration about the effectiveness of the 
proposed feedback-control scheme is whether it is possi- 
ble to select an appropriate set of controller parameters, 
k, that will provide an improvement in the performance 
for a given H p without any a priori knowledge of the 
evolution. We define the increase in success probability 
SP as 



SP 



Pfb - Plin 
Plin ' 



(15) 



speed-up of approximately Ot 



This speed-up occurs 



where Pfb is the success probability for a feedback- 
controlled AQC operation with a given k and total com- 
putation time T and Pn n is the success probability for the 
same AQC operation with linear interpolation over the 
same time T. Figure [4] shows that the average increase in 
success probability due to the feedback-control approach 
as a function of the controller gain is unimodal and that 
there is a clear optimal value of k. These results are aver- 
aged over an ensemble of random Hp's which could repre- 
sent a wide range of computational problems. Therefore, 
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FIG. 4. (Color online) The average percentage increase in 
the success probability (SP) for feedback-controlled interpo- 
lation compared to linear interpolation as a function of the 
controller's gain k. These results are averaged over 10,000 
different 2-qubit problem instances of the form (1131) . 



we can conclude that an appropriate value of fc, which 
will usually provide an increase in P for an arbitrary H p , 
can be found. 



V. CONCLUSIONS 

We have proposed a simple feedback-control scheme 
for solid-state implementations of AQC comprised of flux 
qubits. The scheme makes use of existing on-chip hard- 
ware to perform weak, continuous measurement of the 
ground state curvature. The interpolation speed is dy- 
namically controlled as a function of the ground state 
curvature in order to maximise the success probability. 
The ground state curvature increases to a maximum at 
the minimum ground state gap, where the computation 



speed is reduced to avoid excitation. Then, as the lev- 
els spread back apart and the curvature decreases, the 
computation speed is increased with little risk of excita- 
tion. This approach can be interpreted as attempting to 
optimise the computation speed in a piece-wise fashion 
throughout the entire evolution, compared to calculating 
a single optimal speed based on the limiting minimum 
gap as in 

As a theoretical proof of concept, we have shown that 
a simple proportional gain controller will provide a sig- 
nificant increase in performance when compared to the 
usual linear interpolation scheme for a class of generic 
problem Hamiltonians. Crucially, it was shown that 
this scheme can provide a polynomial speed-up, approx- 
imately 0(n 2 ), in the scaling of the computation time 
with system size, although the exact nature of this scal- 
ing in the limit of large system sizes remains unclear. We 
have also shown that it should be possible to select ap- 
propriate controller parameters that will usually provide 
a significant increase in success probability for an arbi- 
trary problem. We expect that the improvement gained 
from this scheme will be most noticeable in the approx- 
imately adiabatic regime, which is also where we would 
like practical AQC devices to operate, as it offers the 
best compromise between computation time and success 
probability. By tailoring the controller rule, e.g. using 
a more sophisticated non-linear controller rule, or by re- 
stricting ourselves to a more specific class of problems, it 
may be possible to achieve even greater improvements in 
performance using this feedback-control scheme. 
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